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Abstract 

In this paper, we develop a fully discrete Galerkin method for solving initial value fractional 
integro-differential equations(FIDEs). We consider Generalized Jacobi polynomials(GJPs) with 
indexes corresponding to the number of homogeneous initial conditions as natural basis functions for 
the approximate solution. The fractional derivatives are used in the Gaputo sense. The numerical 
solvability of algebraic system obtained from implementation of proposed method for a special case 
of FIDEs is investigated. We also provide a suitable convergence analysis to approximate solutions 
under a more general regularity assumption on the exact solution. 
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1 Introduction 


In this paper, we provide a convergent numerical scheme for solving FIDE 

' X 

V^u{x) = p{x)u{x) + f{x) + A f K(x, t)u{t)dt, x G D = [0,1], 
< 0 

. n(0) = 0, 


( 1 ) 


where q G M"*'P|(0,1). The symbol M'*' is the collection of all positive real numbers. p{x) and f{x) 
are given continuous functions and K{x,t) is a given sufficiently smooth kernel function, u{x) is the 
unknown function. 

Note that the condition n(0) = 0 is not restrictive, due to the fact that ([T]) with nonhomogeneous 
initial condition tt(0) = d, d ^ 0 can be converted to the following homogeneous FIDE 

' _ X 

V'^u{x) = p{x)u{x) + f{x) + X f K{x, t)u{t)dt, x G D = [0,1], 


m = 0 , 


by the simple transformation -u{x) = u(x) — d, where /(x) = /(x) -|- di p{x) + A K{x, t)dt 


I 




Such kind of equations arising in the mathematical modeling of varions physical phenomena, snch as 
heat condnction, materials with memory, combined condnction, convection and radiation problems([3], 

m, m, m)- 

'D‘^u{x) denotes the fractional Caputo differential operator of order q and defines as([8], [l9], [28]) 

P%(x)=XW(x), (2) 

where 

X 


Z^u{x) 



sY ^u{s)ds, 


(3) 


0 

is the fractional integral operator from order ix. T{fx) is the well known Gamma function. The following 
relation holds [8| 


Z^{V^u{x)) = u{x) — u(0). 


(4) 


From the relation above, it is easy to check that © is equivalent with the following weakly singnlar 
Volterra integral eqnation 


u{.) = si.) + X j 


(5) 


Here g{x) = Z^f{x) and K{x, t) = —P(^)+/ — KY, t)ds. From the well known existence 

and uniqueness Theorems([3], [32])) h can be concluded that if the following conditions are fulfilled 




• fix) E c\n), i > 1 

• pix) E c'^(n), i > 1 

• Kix,t) £ C\D), D = {ix,t)]0 <t <x <1}, 1>1 

• K{x, x) Y 0) 


the regularity of the unique solution u(x) of © and also o is described by 

u(x) = ^ + Ui{x; q) E C’'{0, 1] Q C'(n), with |u'(x)| < Cqx'^~^, (6) 

U,k) 

where the coefficients are some constants, Ui{.; q) E (7^(0) and (j, k) := {{j, k) : j,k ^ No, j+kq < 
1}. Here No = N|J{0}, where the symbol N denotes the collection of all natural numbers. Thus, we 
must expect the first derivative of the solution to has a discontinuity at the origin. More precisely, if 
the given functions g{x),p{x) and iF(x,t) are real analytic in their domains then it can be concluded 
that there is a function U = Uizi^Z 2 ) real and analytic at (0,0), so that solutions of ©and also © 
can be written as u{x) = f7(x, x'^)([l], [32]1. 

Recently, several numerical methods for the numerical solution of FIDE’s have been proposed. In 
[25] . fractional differential transform method was developed to solve FIDE’s with nonlocal bound¬ 
ary conditions. In [30] . Rawashdeh studied the numerical solution of EIDE’s by polynomial spline 
functions. In [2], an analytical solution for a class of EIDE’s was proposed. Adomian decomposition 
method to solve nonlinear EIDE’s was proposed in [2T]. In [33], authors solved fractional nonlin¬ 
ear Volterra integro differential equations using the second kind Chebyshev wavelets. In [18], Taylor 
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expansion approach was presented for solving a class of linear FIDE’s including those of Fredholm 
and Volterra types. In [22], authors were solved FIDE’s by adopting Hybrid Collocation method to 
an equivalent integral equation of convolution type. In m. Chebyshev Pseudospectral method was 
implemented to solve linear and nonlinear system of EIDE’s. In [23], authors proposed an analyzed 
spectral Jacobi Collocation method for the numerical solution of general linear FIDE’s. In [TH], au¬ 
thors applied Collocation method to solve the nonlinear EIDE’s. In [2j], Mokhtary and Ghoreishi, 
proved the convergence of Legendre Tau method for the numerical solution of nonlinear FIDE’s. 

Many of the techniques mentioned above or have not proper convergence analysis or if any, very 
restrictive conditions including smoothness of the exact solution are assumed. In this paper we will 
consider non smooth solutions of ([T]). In this case although the discrete Galerkin method can be 
implemented directly but this method leads to very poor numerical results. Thus it is necessary to 
introduce a regularization procedure that allows us to improve the smoothness of the given functions 
and then to approximate the solution with a satisfactory order of convergence. To this end, we propose 
a regularization process which the original equation ([1]) will be changed into a new equation which 
possesses a more regularity properties by taking a suitable coordinate transformation. Our logic in 
choosing proper transformation is based upon the formal asymptotic expansion of the exact solution 
in ([6]). Consider ([T]), using the variable transformation 


X 


VI , V = x'^, t = W‘> , 


w = 


(7) 


we can change ([I]) to the following equation 


M‘^u{v) = p{v)u{v) + 


/» + A 


V 

j K{v, 


w)u{w)dw, 


( 8 ) 


where 


p{v) 


M^u{v) 


1 

p{v'>), f{v) 


1 

fivi), K{v,w) = 


I 




(9) 


From ([6]), the exact solution u{v) can be written as u{v) 

can be easily seen that u'{v) G C(D). It is trivial that for q - 
will be in the form 


--u{v<i) = ^ + Ui{v<i-,q). It 

{pk) 

A, n G N, the unknown function u(v) 


u(v) = ^ + Ui{v^] q), n G N, 

ihk) 


which is infinitely smooth. Then we can deduce that the solution u{v) of the new equation ([8]) 
possesses better regularity and discrete Galerkin theory can be applied conveniently to obtain high 
order accuracy. 

In the sequel, we introduce the discrete Galerkin solution un{v) based upon GJPs to ([8]). Since 
the exact solutions of ([I]) can be written as u{x) = u{v) then we define un^x) = U]\fiv), x,v G Q as 
the approximate solution of ([!]). 

Spectral Galerkin method is one of the weighted residual methods(WRM), in which approximations 
are defined in terms of truncated series expansions, such that residual which should be exactly equal 
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to zero, is forced to be zero only in an approximate sense. It is well known that, in this method, 
the expansion functions must satisfy in the supplementary conditions. The two main characteristics 
behind the approach are that, first it reduces the given problems to those of solving a system of 
algebraic equations, and in general converges exponentially and almost always supplies the most terse 
representation of a smooth solution([6], [E], [HT]). 

In this article, we use shifted GJPs on Q, which are mutually orthogonal with respect to the shifted 
weight function = (2 — 2u)“(2u)^ on II where a, 13 belong to one of the following index sets 

Ml = {(a, /3); a,/? < —1, a,l3& Z}, A 2 = {(a,/3);a < —1,(3 > —1, a & Z, (3 & M}, 

A /3 = {{a,(3)]a > —1,(3 < —1, aEM,/3GZ}, A /4 = {(a,/3);a,/3 > —1, q:,/3gM}, 

where the symbol Z is the collection of all integer numbers. The main advantage of GJPs is that these 
polynomials, with indexes corresponding to the number of homogeneous initial conditions in a given 
FIDE, are the natural basis functions to the Galerkin approximation of this problem/ 112]. |13jl. 

The organization of this paper is as follows: we begin by reviewing some preliminaries which are 
required for establishing our results in Section 2. In Section 3, we introduce the discrete Galerkin 
method based on the GJPs and its application to ([8|). Numerical solvability of the algebraic system 
obtained from discrete Galerkin discretization of a special case of ([8]) with 0 < (? < ^ and p{v) = 1 
based on GJPs is given in Section 4. Convergence analysis of the proposed scheme is provided in 
Section 5. Numerical experiments are carried out in Section 6. 


2 Preliminaries and Notations 


In this section, we review the basic definitions and properties that are required in the sequel. 
Defining weighted inner product 


Ui,U2 

and discrete Jacobi-Gauss inner product 


q,/3 


f ui{v)u2{v) 6'^’^{v)dv, 

Jn 


\ ^ 

ui,u2] 

J N,a,P 


we recall the following norms over D 


l^lla,/3 = \U,U 


a,l3 


l7/l|2 


g=(u,u) , ||u||oo = SUp|u(u)|. 


Here, and are the shifted Jacobi Gauss quadrature nodal points on D and corresponding 
weights respectively. 

The non-uniformly Jacobi-weighted Sobolev space denotes by ^{Q) and defines as follows 

< 00 ; 0 < s < A:}, 

equipped with the norm and semi-norm 

k 

I I'^l \a,0,k — / J P lla-|-s„S+s> \u\a,/3,k = \\U \\a+k,/3+k- 

s=0 
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The space ^{Q) distinguishes itself from the usual weighted Sobolev space ^{Q) by involving 
different weight functions for derivatives of different orders. The usual weighted Sobolev space 
is defined as 


= {u : < oo; 0 < s < A;}, 


equipped with the norm 



k 


E 




||2 


s=0 


We denote the shifted GJPs on Q. by Gn'^{v) and define as 

(2 - {a, 13) e Ml, no = -{a + 13), 

(2 - 2'(;)-"J"";^('(;), (a,/3) E M 2 , m = -a. 


{2v) ^Jn'-hv), 


7"’^ iui 

V ^n—nQ\^)'> 


{a,(3)eM3, no = -f3, 

{a,(3)GM4, no = 0, 


( 10 ) 


where Jn'^{v) is the classical shifted Jacobi polynomials on 11; see m- An important fact is that 
the shifted GJPs {Gn'^{v);n > 1} form a complete orthogonal system in see([l2], [13]). 

To present a Galerkin solution for ([8|) it is fundamental that the basis functions in the approximate 
solution satisfy in the homogeneous initial condition. To this end, since Gn~^(0) = 0, n > 1, then 
we can consider {Gn~^(u), n > 1} as suitable basis functions to the Galerkin solution of Q. 

From (IlOp and the following formula [9] 


ja,/ 3 f , f r(f +/3 + l)r(i + A: + a + /3 + 1) l. 0 ^ 

j, T{k + l3 + l)ni + a + l3 + l){i-ky.k< ” ■ 


A:=0 


we can obtain the following explicit formula for ^(u) 


2—1 


= (2^)JS(!.) = 2 ^^(-1)*-' 


-l-k 


{i + k)\ 


k+l 


k=0 


{k + l)\{i - 1 - k)\k\ 


i > 1. 


( 11 ) 


For any continuous function Z{v) on H, we define the Legendre Gauss interpolation operator Zn, 


as 


N 


InZ{v) = 




s=0 


I 70>0||2 


A7Mj0,0(^)_ 


( 12 ) 


TV, 0,0 


Let Vn be the space of all algebraic polynomials of degree up to N. We introduce Legendre 
projection Utv : M{kl) —)• Vn which is a mapping such that for any Z{v) E M{Q), 


Z — UnZ, 4> ] =0, \/(j) E Vn- 

0,0 


(13) 
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3 Discrete Galerkin Approach 


In this section, we present the numerical solution of ([8]) by using the discrete Galerkin method based 
on GJPs. 

Let 

N 

(14) 

i=l 

be the Galerkin solution of ([8]). It is trivial that iiAr(O) = 0. 

Galerkin formulation of ([8|) is to find UNiv), such that 

(M‘^UN,G^rA = (piv)uN{v),G^r^] +(f{v),G^p-A +x(lC{uN),G^p-^] , i = l,2,...,N 

V /o,-i V /o,-i V / 0,-1 V /o,-i 

(15) 

V 

where JC{un) = f K{v,w)uN{w)dw. 

0 

Applying transformation w{0) = v6, 0 G hi we get 

1 

IC{un) = fCe{uN) = V j K{v,w{9))uNiw{9))d9. (16) 

0 


Substituting (fT6l) in (fTSl) yields 


0,-1 


M'‘uN,Gf M = {p{v)uN{v),Gf +{f{v),G' 


0,-1 


0,-1 


+ x(lCeiuN),G^p-^ 

0,-1 V / 0,-1 

i = l,2,...,N. (17) 


Inserting (fH)) in (fT7)l we get 


= , i = l,2,...,N. (18) 


0,-1 


Following the relation ^(i;)5'^’ ^(u) = G°g\(i;), we can rewrite m as 


= (/(!>),G"i‘i) , i=l,2,...,JV. (19) 

0,0 
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Now, we try to find an explicit form for A4'^G^’ To this end, using dill) we have 






T(l-q) 

2 


k=0 
i-1 


-1-k 


(j + k)l 


(k + - 1 - k)l 


(vi—wi) (20) 








{yi — wi) '^w^dw. 


Applying the relation [TT] 

f f - ki (qTTCSc{'Kq)T{q +qk)\ ^ „ 


in d20l) we can obtain the following explicit formula for M^G' 


0 ,- 1 . 


9 i-i 

M^^G^-^v) = (-l)J 


r(l-g) 


_l_fc {j + k)\ (qT:csc{'Kq)T{q + qk)\^k ( 2 I) 


k=0 


{k\f{j - l-k)\\ T{q)T{l + kq) 


Substituting (pTIl in (fTOll we obtain 


E<',j('I'MW.G?ti) -(pWG“-‘W,G«i\) -A(c„(G”'-‘),Gyi) I 

j=l IV /o,o V /o,o V /o,oJ 

= (f{v),G°Y) > i = l,2,...,iV. (22) 


0,0 


In this position, we approximate the integral terms of (|22l) using (N + 1)—point Legendre Gauss 
quadrature formula. Our discrete Galerkin method is to seek 


N 


Un{v) = 


(23) 


2 = 1 


such that coefficients satisfies in the following algebraic system of linear equations 

E“l('f«('’).G"i‘i) -(pWgJ'-'w.gE) -a(k„.„(g“'-‘),g“i‘,) I 

j=l IV /0,0 V / Nfi,0 V /Af,0,oJ 

= = (24) 

where 

N 

ICnAG^Y) = vYKi'>J^H^k))GY\wi9k))6k 


N, 0,0 


(25) 


k=0 


Here shifted Legendre Gauss quadrature points on H and correspond¬ 

ing weights respectively. Note that, from ()2ip we can see that 'l'j^g(u) is a polynomial from degree at 

most N, then we have ( 'l'j^g(u), j = ( ^j^q{v), j . It is trivial that the solution of (IMl) 
gives us unknown coefficients in 
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4 Existence and Uniqueness Theorem for Discrete Galerkin Alge¬ 
braic System 

The main object of this section is providing an existence and uniqueness Theorem for a special case 
of the discrete Galerkin algebraic system of equations ([Ml) with p{v) = 1 and 0 < q < ^. Throughout 
the paper, Ci will denote a generic positive constant that is independent on N. 

First, we give some preliminaries which will be used in the sequel. 

Definition 4.1 Let X^y be normed spaces. A linear operator A : X ^ y is compact if the set 
{Ax\ \\x\\x < 1} has compact closure in y. 


Theorem 4.2 \1(A 11, "3^ Assume X, y he two banach spaces. Let 


V = Av + f, (26) 

be a linear operator equation where A : X ^ y is a linear continuous operator, and the operator I — A 
is continuously invertible. As an approximation solution of ^26\) we consider the equation 


VN = AnVn + ^Nf, 


(27) 


which can be rewritten as 

VN = BnAvn + Snvn + ^Nf, (28) 

where An is a linear continuous operator in a closed subspace y of y. Bn, Bn : T —>■ T are linear 
continuous projection operators and Sn = An — Bn A is a linear operator in y. If the following 
conditions are fulfilled 


• for any Z G y we have ||5Ar.^|| —>-0 as N ^ oo 

• ||.4 — .SAr.4|| —>• 0 as N ^ oo 

• \\f - BNfW ^ 0 as N ^ oo 

then Hgj) possesses a uniquely solution vn £ T, for a sufficiently large N. 


Lemma 4.3 Let X,y be banach spaces and y be a subspace ofy. Let Bn : y ^ y be a family of 
linear continuous projection operators with 

BnU y as N ^ oo, y G y. 


Assume that linear operator A : X ^ y be compact. Then 

||.4 — .SAr-4.|| —0 as N ^ oo. 

Lemma 4.4 (Interpolation error bound, \31^ ) LetXNZ be the interpolation polynomial approximation 
of the function Z{v) defined in For any Z{v) G Bqq{Q) with k > 1, we have 


Z-lNZ\\o,o<CN-’^\Z\o^o^k- 




Lemma 4.5 For every bounded function Z{v), there exists a constant C independent of Z such 
that 


sup ||Xjv.Z’||o,o < C sup \Z{v)\. 

N ’ V 

Lemma 4.6 (Legendre Gauss quadrature error bound lSl^ ) If Z{v) G Bqq{LI) for some k > 1 and 
$ G Vn, then for the Legendre Gauss integration we have 


(Z, ‘^)o,o - {Z, 4>)Ar,o,o 


< CN ^||Z||o,o,A:||‘^’||o, 0 - 


Now we intend to prove existence and uniqueness Theorem for a special case of the discrete Galerkin 
system (|Mll with p{v) = 1 and 0 < g < ^. 

Theorem 4.7 (Existence and Uniqueness)Let 0 < q < ^ and p{v) = 1. // ^ has a uniquely solution 
u{v) then the linear discrete Galerkin system [2^ has a uniquely solution UAr(u) G Vn for sufficiently 
large N. 

Proof: Our strategy in proof is based on two steps. First, we try to represent (I24h in the operator form 
(I28p . Then by applying Theorem 14.21 to operator form obtained in the first step the desired result 
have been concluded. 

Step 1: In this step, we show that the discrete Galerkin system (j24l) can be written in the operator 
form (f28]l . To this end, consider ([8|) and define 

flNiv) = M'^UNiv) - UNiv) - XICN,e{uN) - f{v). 


According to the proposed method, we have 

nN{v),G^p\{v)] =0, i = l,2,...,iV. 

/ AT,0,0 

From interpolation and Legendre Gauss quadrature properties, we can write 


=[ZN{nN),G(f^{v)] =[ZN{nN),G(f^{v)] =0, i = l,2,...,iv. 

/ Af,0,0 V / AT,0,0 V / 0,0 

(29) 

Since Zn{1Zn{v)) is a polynomial, it can be represented by a linear orthogonal polynomial expan¬ 
sion based on as 


■v0,l 


■< 0,1 


Zn{Vn) = ^ 


^ \ZN{nN),G^i 


0,-1 




2=1 


\G) 


0,-l||2 




0,0 ^0,-1 


110,-1 .=1 l|Gr 

Using the relations (f29]l and (pOl) yields Zn{1Ln) = 0. Thus 


'0,-l||2 
lo,-i 


Gr\v). (30) 


Zn { M^un — un{v) — XKn,9{un) — f ] = 0) 
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which can be rewritten as 


un{v) = In 




Un 


^nI 


XnTnUn — Xj^f , 


and thereby 


un{v) = Tii^TuN + SnUn - X^f, 


where X^r and IIjv are defined in (fl^ and (fOl) respectively and 


r = 

Tn = — XKn^q, 

Sn = XnTn - HatT. 


(31) 

(32) 


Since ()32p is the form in which the discrete Galerkin method is implemented, as it leads directly 
to the equivalent linear system (|24p . It can be easily check that (|32ll can be considered in the operator 
form (I28p by assuming 

A' = y = L\n), y = VN, 

vn = Un, Bn = ^n, BN=X]y, (33) 

A = T, An = XnTn, 

which can be completed the desired result of step 1. 

Step 2: In this step we intend to apply Theorem 14.21 with the assumptions (|33p to prove the 
Theorem. To this end, following Theorem 14.21 we must show that 

1) for any Z E Vn we have ||5ArZ||o,o = \\XnTnZ — nArTX||o,o —)• 0 as iV —)• oo, 

2) ||7~ — nTvT^llojO —^ 0 as —y oo, (34) 

3) 11 / ~ ^N / ||o,o 0 as N —7> oo. 


First, we prove the first condition in (I34p . For this, we can write 


|‘5Ar^||o,o = II (^nTn — ■Z’llo.o < II^^Af (tn — XHo^o + || (xn — nAr^T.Z'||o,o- (35) 


Since M.‘>Z € Vn for any Z E Vn then 


^TvXllo.o < A IIXatI )CN,e — )-2^110,0 + II { — Htv )/Ce)Z||o,o 


(36) 


Using Lemmas 14.51 and 14.61 and the relations (jl6p and (1251) we can obtain 

|| 2 lArf/CAr,(?-/C 0 ^z||o,o < sup|/C 0 (z(u)) - ICn, 9 {z{v))\ < CiN~^^ sup r||iF(n, u;( 6 »))||o,o,fci lb^(w'(6'))||o,oY 
V / v£fl vGfi V / 

(37) 

where norms ||7L(u, t(;( 0 ))||o,o,A:i and ||nZ(r(;(6())||o,o are applied with respect to the variable 0. 
According to Lemma 14.61 we have 


0,0 V / TV,0,0 


< CN->^^\\}Ce\\o,oM\\Js’X,o, s = 0,1,...,fV, 


(38) 
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where norm ||/C 6 )||o,o,fe 2 i® applied with respect to the variable v. Using ([381) and the relations (fT^ . 
(USD we can yields 

(39) 


{^N ~ ^ 0 as N —> oo. 


Substituting ([37)1 and ([3^ in (f36]l we can conclude the first condition in ([3^ . Applying Lemma 
14.61 gives us the third condition in ([34|1 . To complete the proof, it is sufficient that we prove the second 
condition of ([M|l . To this end, we apply Lemmawith the assumptions (l33l) . Since ||y — nAry||o,o —t 
0 as A" ^ oo for y G L^(n)(see[]), the second condition in (IMll can be achieved by proving compactness 
of the operator T. Since k{v, w) is a continuous kernel, then the operator T will be compact if : 

Aq o(f^) —^ L^(n) be a compact operator. For this, dehne M = [/q^ ~ ) dwdv] ^ < oo. 


For u G FrQQ(n), using Cauchy-Schwartz inequality we have 


ll-^'^^llo,o = J u'{w)dw\'^dv < j 


< 


ni: 


yq — WQ 


-Q 


u'{w)\dw ) dv < 


< M^\\u'\\lo< M^\\u"^ 


Koi^y 


1 1 
yq — yjq 


-Q \ 2 

u\w)\dw ) dv 


1 i 

yq — W 1 


-2q 


dw 


w 


{w)\^dw^dv 


So M.^u G L^(n) and ||A4'?|| = sup <1 


u 


7 ^ 0,u G j < M < oo. Therefore 

A4'^ : Hqq{Q) —)• L^(0) defined by ([9]), is a bounded operator. If we proceed same as Theorem 3.4 in 
m, we can conclude compactness of the operator A4'^(u). Then the operator T is a compact operator 
and from Lemma 14.31 we have 

ll'T — Hjy'TIIq^q —y 0 as A —y oo 

In this position, all conditions in (IM[) that are required to deduce existence and uniqueness of 
solutions of the discrete Galerkin system (1241) have been proved and then the proof is completed. □ 


5 Convergence analysis 

In this section, we will try to provide a reliable convergence analysis which theoretically justify con¬ 
vergence of the proposed discrete Galerkin method for the numerical solution of a special case of ([ 8 ]) 
with p{v) = 1 . 

In the sequel, our discussion is based on these Lemmas: 

Lemma 5.1 jdl^ For any Z G Hqq)!!), we have 

\\InZ\\o,o<c(\\Z\\o,o + N-^\\Z'\\i^X 


Lemma 5.2 The fractional integral operator is hounded in and 


||X^Z||o,o < CIlZlIo.o, ^ G L\n). 
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Lemma 5.3 ^(Gronwall inequality)Assume that Z{v) is a non-negative, locally integrable function 
defined on ri which satisfies 


Z(v) < b{v)-\-B / [v — Z{w)dw, w G 

Jo 

where a,fi> —1, b{v) > 0 and B >0. Then, there exist a constant C such that 

pv 

Z(v) < b{v)C / {v — wPw^b{w)dw, wGQ. 


Now, we state and prove the main result of this section regarding the error analysis of the proposed 
method for the numerical solution of a special case of ([ 8 ]) with p{v) = 1 . 

Theorem 5.4 (Convergence)Let u{x) and u{v) are the exact solutions of the equations ^f) and ^ 
respectively that is related by u{x) = u{v). Assume that un{x) = UNiv) be the approximate solution 
of (QP, where un{v) is the discrete Galerkin solution of the transformed equation If the following 
conditions are fulfilled 

1. f{v) E -Bg o(n) for ki>l, 

2. /C(tt) E Bq)q{H) for k 2 > 1, 
then for sufficiently large N we have 

||eAf(^i)||o,o < <^-^”^^1/10,0,^! + l^(^)lo,o,fc2^ 
where = min{A:i, ^ 2 } and eM{u) = u{x) — un{x) denotes the error function. 


Proof: 


Since AA‘^un is a polynomial from degree of at most N then we can rewrite (131 p as 


— Un{v) —Tn( ] =X7v/. 


(40) 


Subtracting ([ 8 ]) from (1401) and some simple manipulations we can obtain 


M'^cm = eN{v) + exj,rif) + A/C(eAr) + )^ex^{LJ{uN)) + ZlN\^e{uN) - l^Nfi{uN)j , (41) 

where 6 x^( 5 ) = g —X]x{g) denotes the interpolating error function and eM{v) = u{v) — um{v) is the 
error function of approximating ([ 8 ]) using discrete Galerkin solution un{v). 

Applying the transformation ([7]) to the operator we get 


M\u = 


m 


1 1 

VI — wi 


q-l --1 

u{w) - dw. 


Following dH) it is easy to check that 

A4f (A4 ^m(u)) = u{v) — u{0). 


(42) 
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Applying the operator on the both sides of (HIT) and using (|12]) we can yield 


\eN{v)\ <-A4i(|eAr| + |A/C(eAr)|) + 


■Ml ^exffif) + Xexj^{lC{uN)) + XIn (^e{uN) — ,e{uN) 


■ (43) 


Since 


^ L/ Ciy ^ ^ 

Mf J \K{v,w)e]\f{w)\dw = ^ J J (y'i — w'^Y -|^(tc, s)||eAr(s)|dsdt(; 

0 0 0 


m 


[vi — Wl) 






-\K{w, s)||eAr(s)|(ir(;fis 


0 s 


r(g) 


V 

J Ki{v,s)\eN{s)\ds < 


l|-^l(^.g)l|c 

r(g) 


|eAr(s)|ds, (44) 


then (j33]) can be rewritten as 


V 1_^ V 

\eN{v)\ < ^ \eNiw)\dw + |A|||^i(n, s)||oo J \eNiw)\d- 

0 0 


w 


+ 


■Ml[exjy(f) + XeXff{1C{uN)) + XlN\K,e{uN) - K^N,e{uN) 




V 

< — w'^Y ^|eAr('W^)|d'(i'^ 


(45) 


+ 


■Mf ( exj^if) + Xexj^{JC{uN)) + XlN^JCeiuN) - ICN,e{uN) 


where 


/* / 1 lx n—^ IC ~ ~ f 

Ki{v,s) = J [v'l—w’iY - \K{w,s)\dw, (7 = max I 


1 |A|||Ki(n,s)||, 


9r(9)’ T{q) 


||l+(n'J —to'j) 




Using Gronwall ineaualitvlLemma 15.3p in (I45p yields 


ewllo 1-1 < Cl I ||A4^ I exj^if) + Aei^(/C(MAr)) + AXatI ^^^{un) — ) ) |lo i-i ) • (46) 


It can be checked that by applying the transformation ([6]) we obtain 


|eiv(t^)|lo i_i = g22 ie7v(x)||o,o- 


(47) 


On the other hand, by applying Q and Lemma 15.21 we can get 


\\MlZiv)\\l,_^ = q2-^\\I^Zix)\\l, < Cq2^-YZ{x)\\l, = C||Z(i;)||^_,_^ 


--ii 


(48) 
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where Z{x) is a given function and Z{v) = z{v^). Using the relations (H7)) and (|i8]) in (H6]) we can 
obtain 


q 2 <i ^||e7v||o,o < C 2 \\exj^{f) + \ext^{lC{uN)) + XIiq\^e{uN) — ^N,e{uN)j\\Q^i_i 

< <^2 ^||e2;jv(/)llo,i-i + l■^ll|exJv(^('“^^))llo,|-l + \M\\^n (^}C0{un) - JCisf^eiuN)^ ||o,i 
Since then we have ||.||q i_^ < ||.||o,o and thereby ()39]l can be written as 

||eAr||o ,0 < ^3 (/il +/i 2 + ^3 ) , (50) 


where 


= ||exjv(/)l|o,o, 1^2 = ||exjv(^(^Af))l|o,o, -^3 = \\ZN{^e{uN) - I^N,e{uN))\\o,o- 


Now, it is sufficient to find suitable upper bounds to the above norms. According to Lemma 14.41 
we have 


(51) 


= ||exjv(/)llo,o < (^3-^ M/lo,o,fci, 

-^2 = ||exjv(^('“Af))llo,o < C'4A^~^^|Aii(u — eAr(u))|o,o,fc2' 

Using Lemmas 14.51 and 14.61 we can obtain 

^3 = W^Nil^eiuN) — )(^N,e{uN)) ||o ,0 < sup \ICg{uN) — ICN,e{uN)\ 

vGQ. 

< C7N~’^^sup r||.A(u,u;(6»))||o,o,fc3ll^^Ar(w'(6'))||o,o 

vefi V 

< C7N~^^f sup||.A(u,u;(6»))||o,o,fc3) ll'«Ar(w')||o,o 

V veQ / 

< C'7Ai"*^3^sup||.A(u,u;(6»))||o,o,fc3^ (^ll^llo,o + ||eAr||o,o^(52) 

where norms \\K{v,w{ 6 ))\\ofi^k 3 and ||unAr(t(;(0))||o,o are applied with respect to the variable 9. 

For sufficiently large N the desired result can be concluded by inserting (1511) and (f52]l in (f50l) . □ 


6 Numerical Results 

In this section we apply a program written in Mathematica to a numerical example to demonstrate the 
accuracy of the proposed method and effectiveness of applying GJPs. the ’’Numerical Error” always 
refers to the L^-norm of the obtained error function. 

Example 6.1 Consider the following FIDE 

X 

T>^u{x) = u{x) + f {x) + ^ j y/xiu{t)dt, u(0) = 0, 

0 
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with the exact solution u(x) = and 

^ ’ \/X 




fix) = 


—X + cos X + ) cos I — Ji(f) sin | ) | — 2sinx 

2-v/x 


/Tere Jn(^) gives the Bessel function of the first kind. 


This example has a singularity at the origin, i.e., 


u' (x) ~ 


1 



In the theory presented in the previous section, our main concern is the regularity of the trans¬ 
formed solution. For the present problem by applying coordinate transformation x = the infinitely 
smooth solution 


u{v) = u{v‘^) = 


smu 


V 

is obtained. The main purpose is to check the convergence behavior of the numerical solutions with 
respect to the approximation degree N. Numerical results obtained are given in the Table 1 and Figure 
1. As expected, the errors show an exponential decay, since in this semi-log representation (Figure 1) 
one observes that the error variations are essentially linear versus the degrees of approximation. 


Table 1: The numerical errors of example O 


N 

Numerical Errors 

2 

7.86 X 10"^ 

4 

4.71 X 10"^ 

6 

2.15 X 10-6 

8 

1.05 X 10-® 

10 

4.42 X 10-1° 

12 

2.05 X 10-^2 

14 

2.64 X 10-^“^ 

16 

1.91 X 10-16 
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N 


Figure 1: The numerical errors of example 16.11 with various values of N. 
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